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METHOD AND APPARATUS FOR FREQUENCY OFFSET ESTIMATION, 
AND SYSTEM UTILIZING SAME 



Field of the Invention 

The present invention relates to frequency offset estimation, and 
more particularly to frequency offset estimation in the presence of direct 
current offset. 

Background of the Invention 

With the advent of personal wireless communication devices, there 
is a need for small portable units that have low power dissipation, and 
can be provided at relatively low cost. To meet these requirements, there 
is interest in a transceiver with direct conversion architecture (DCA). 
Different from conventional heterodyne receivers, the receiver of a DCA 
recovers base band information directly from the carrier signal without 
passing through an intermediate frequency stage and image signal 
rejection. This is accomplished by translating the radio frequency signal 
directly to zero frequency, and then employing a low pass filter to 
suppress any high frequency interference. The DCA receiver has several 
advantages over the heterodyne receiver, and meet the requirements of 
small form factor, low power dissipation, and low cost. However, a 
disadvantage is the down converted base band extends to zero 
frequency, and extraneous offset voltages, referred to as direct current 
(DC) offset, can corrupt the base band signal, and saturate subsequent 
base band processing stages. In addition, a DCA receiver also suffers 
I/Q mismatch, a result of mismatches between amplitude of the I and Q 
signal during quadrature down conversion. There are also concerns 
with even-order distortion, flicker noise and local oscillator leakage that 
adversely affect the performance of a DCA receiver. 



Orthogonal frequency division multiplexing (OFDM) is a bandwidth 
efficient modulation technique that is used for high-speed wireless data 
transmissions. OFDM transmits data on densely packed orthogonal 
sub-carriers that are spaced in frequency exactly at the reciprocal of 
symbol interval. A primary advantage of OFDM is it is utilized for 
transmission at relatively low complexity in multipath channels. A 
disadvantage of OFDM systems is its sensitivity to frequency offset. 
Frequency offset can cause a reduction of signal amplitude at the 
output of OFDM demodulator, and also introduces inter-carrier 
interference (ICI) from the other sub-carriers, thus destroying the 
orthogonality of OFDM. 

Combining the DCA and OFDM provides a transceiver that has a 
small form factor, has low power dissipation, is low cost and is 
bandwidth efficient. However, the DC offset and the frequency offset 
must be addressed. Typically, the approach that is taken is to estimate 
the DC and frequency offsets and then apply compensation based on 
the estimated offsets. In data transmission, offset estimation is typically 
performed during the reception of the preamble of each data packet, 
because in a packet switched system each data packet can be treated 
independently. 

With reference to FIG. 1 a preamble 100 is shown with a time line 
t, where the preamble reception starts at tO, and ends at t3. During a 
first time interval tO to tl a first portion 105 of the preamble is used for 
DC offset estimation; during the second time interval tl to t2 a second 
portion 110 of the preamble is used for DC offset compensation; and 
during a third time interval t2 to t3 a third portion 115 is used for 
frequency offset estimation. Here we assume that the receiver amplifier 
is not saturated by the presence of the DC offset. This conventional 
method sequentially estimates the DC offset, compensates for the DC 
offset, and then estimates the frequency offset. Of course it assumes 
that the DC offset estimation method is based on the initial reception of 



the first portions 105 of the preamble 100, and that the information 
from the first portion 105 is available. 

Due to the limited time available to process the preamble 100, tO to 
tl interval for DC offset estimation, and tl to t2 interval for DC offset 
compensation, one disadvantage is that residual DC offset may still be 
present when the frequency offset is later estimated in the time interval 
t2 to t3. Consequently, the residual DC may result in inaccurate 
frequency offset estimation. In addition, since part of the reception time 
of the preamble 100, interval tO to t2, is used for DC offset estimation 
and compensation, the available time for frequency offset estimation, 
the interval t2 to t3, is reduced. Hence, a further disadvantage of this 
method is the short time during which frequency offset estimation is 
performed, and the consequent inaccuracy of the frequency offset 
estimation. 

Brief Summary of the Invention 

The present invention seeks to provide a method and apparatus for 
offset estimation and system utilizing same, which overcomes or at least 
reduces the abovementioned problems of the prior art. 

Accordingly, in one aspect, the present invention provides a data 
communication receiver for receiving at least one data packet, wherein 
the at least one data packet has a variety of inherent offsets, the data 
communication receiver comprising: 

an input coupled to receive the at least one data packet; and 
a plurality of offset estimators coupled to the input, the plurality 
of offset estimators comprising: 

at least one of the plurality of offset estimators for 
estimating at least one of the variety of inherent offsets from at least a 
received portion of the at least one data packet while the at least one 
data packet is being received, and the at least one of the plurality of 



offset estimators having an output for providing an estimate of the at 
least one of the variety of inherent offsets; and 

at least another one of the plurality of offset estimators 
coupled to the output of the at least one of the plurality of offset 
estimators, the at least another one of the plurality of offset estimators 
for estimating at least another one of the variety of inherent offsets from 
the at least the received portion of the at least one data packet while the 
at least one data packet is being received and from, the estimate of the 
at least one of the variety of inherent offsets, and the at least another 
one of the plurality of offset estimators having an output for providing 
an estimate of the at least another one of the variety of inherent offsets. 

In another aspect the present invention provides a data 
communication receiver for receiving data packets, the data 
communication receiver comprising: 

an input coupled to receive at least one data packet, wherein the 
at least one data packet has a variety of inherent offsets; 

at least one offset estimator for estimating at least one of the 
variety of inherent offsets from at least a received portion of the at least 
one data packet while the at least one data packet is being received, 
and the at least one offset estimator having an output for providing an 
estimate of the at least one of the variety of inherent offsets; 

at least one compensator coupled to the input for receiving the at 
least one data packet, and coupled to the output of the at least one 
offset estimator for receiving the estimate of the at least one of the 
variety of inherent offsets, the at least one compensator for 
compensating the at least one of the variety of inherent offsets of the 
data packet using the estimate of the at least one of the variety of 
inherent offsets to produce at least a portion of a partially compensated 
data packet, and wherein the at least one compensator has an output 
for providing the at least the portion of the partially compensated data 
packet corresponding to the at least one data packet; 
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a plurality of offset estimators coupled to the output of the at 
least one compensator, the plurality of offset estimators comprising: 

at least one of the plurality of offset estimators for 
estimating the at least one of the variety of inherent offsets from the at 
5 least the portion of the partially compensated data packet while the at 
least the portion of the partially compensated data packet is being 
received, and the at least one of the plurality of offset estimators having 
an output for providing an estimate of the at least one of the variety of 
inherent offsets; and 
10 at least another one of the plurality of offset estimators 

coupled to the output of the at least one of the plurality of offset 
estimators, the at least another one of the plurality of offset estimators 
Nj for estimating at least another one of the variety of inherent offsets from 

fS the at least the portion of the partially compensated data packet while 

15 the at least the portion of the partially compensated data packet is 
being received, and from the estimate of the at least one of the variety of 
inherent offsets, and the at least another one of the plurality of offset 
estimators having an output for providing an estimate of the at least 
another one of the variety of inherent offsets, 
20 In yet another aspect the present invention provides a method for 

offset compensation comprising the steps of: 

a) determining at least one data packet is being received; 

b) estimating direct current offset of the at least one data packet to 
produce an estimate of direct current offset; 

25 c) estimating frequency offset of the at least one data packet using 
the estimate of direct current offset; 

d) compensating the at least one data packet using the estimate of 
direct current offset and the estimate of frequency offset to produce a 
compensated data packet; 
30 e) providing the compensated data packet corresponding to the at 
least one data packet. 
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In still another aspect the present invention provides a method for 
frequency offset estimation comprising the steps of: 

a) determining samples of at least one data packet are being 
received; 

b) determining direct current offset power using the samples of the 
at least one data packet; 

c) adding a predetermined constant to the direct current offset 
power; 

d) storing a predetermined number of the samples of the at least one 
data packet; 

e) performing complex conjugation on the predetermined number of 
the samples to produce a result; 

f) multiply the result of the complex conjugation with the samples 
of the at least one data packet to produce an output; 

g) determine the average of the output from step (f) by another 
predetermined number to provide an average; 

h) determine argument of the average and divide the argument by 
the negative value of the predetermined number to produce a resulting 
quotient; and 

i) multiply the resulting quotient by the sum of the constant and 
the direct current offset power determined in step (c) to produce a 
frequency offset estimate. 

In yet still another aspect the present invention provides a 
method for frequency offset estimation comprising the steps of: 

a) determining samples of at least one data packet are being 
received; 

b) determining direct current offset power using the samples of the 
at least one data packet; 

c) storing a predetermined number of the samples of the at least one 
data packet; 

d) performing complex conjugation on the predetermined number of 
the samples to produce a result; 
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5 Brief Description of the Drawings 



e) multiply the result of the complex conjugation with the samples 
of the at least one data packet to produce an output; 

f) determine the average of the output from step (e) by another 
predetermined number to provide an average; 

g) combine the direct current offset power from step (b) with the 
average from step (f) by subtraction; 

h) determine argument of the result of subtraction from step (g) and 
divide the argument by the negative value of the predetermined number 
to produce a frequency offset estimate. 
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An embodiment of the present invention will now be more fully 
described, by way of example, with reference to the drawings of which: 



H " 15 FIG. 1 shows a time-activity diagram of a prior art offset estimation 

^ and compensation system; 

Hi FIG. 2 shows a functional block diagram of one embodiment of an 

:J; offset estimation and compensation system in accordance with the 

fU present invention; 

20 FIG. 3A shows a time-activity diagram of the offset estimation and 

compensation system in FIG. 2; 

FIG. 3B shows an alternate time-activity diagram of the offset 
estimation and compensation system in FIG. 2; 

FIG. 4 is a flowchart detailing the operation of the offset estimation 
25 and compensation system in FIG. 2; 

FIG. 5 shows a functional block diagram of another embodiment of 
an offset estimation and compensation system in accordance with the 
present invention; 

FIG. 6A shows a time-activity diagram of the offset estimation and 
30 compensation system in FIG. 5; 

FIG. 6B shows an alternate time-activity diagram of the offset 
estimation and compensation system in FIG. 5; 
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FIG. 7 is a flowchart detailing the operation of the offset estimation 
and compensation system in FIG. 5; 

FIG. 8 shows the structure of a preamble of a data packet; 
FIG. 9 is a model of a communication channel on which data 
5 packets having the preamble in FIG. 8 are transmitted; 

FIG. 10 provides a graph showing the deterioration in frequency 
offset estimates with increasing DC offset power; 

FIG. 11 shows a functional block diagram of a prior art frequency 
offset estimator without averaging; 
10 FIG. 12 shows a functional block diagram of another prior art 

frequency offset estimator when averaging is employed; 
P FIG. 13 shows a functional block diagram of a first embodiment of 

p a frequency offset estimator in accordance with the present invention; 

N FIG. 14 is a flowchart detailing the operation of the frequency offset 

m 15 estimator in FIG. 13; 
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FIG. 15 shows a functional block diagram of a second embodiment 
of a frequency offset estimator in accordance with the present invention; 
and 

FIG. 16 is a flowchart detailing the operation of the frequency offset 
20 estimator in FIG. 15; 

FIGS. 17A-F provide graphs showing the performance of the first 
and second embodiments of the present invention based on the 
preamble of a smaller number of data packets; 

FIGS. 18A-F provide graphs showing the performance of the first 
25 and second embodiments of the present invention based on the 
preamble of a larger number of data packets; and 

FIGS. 19A-C provide graphs comparing the performance of the 
second embodiment of the present invention to a prior art frequency 
offset estimator. 

30 
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Detail Description of the Drawings 

A direct current (DC) and frequency offset estimation and 
compensation system receives data packets, each having a preamble. 
5 The system comprises a DC offset estimator and a frequency offset 
estimator that operate concurrently to produce estimates of DC and 
frequency offsets of the data packet, which are determined from a 
portion of the preamble as it is received. As the rest of the preamble is 
received, or alternatively after the preamble is received, a compensator 

10 receives the estimates of DC and frequency offsets and compensates the 
remaining portion of the data packet to produce a DC and frequency 
offset compensated data packet. 

The concurrent operation of the DC and the frequency offset 
compensators advantageously allows more time to produce the 

15 estimates of DC and frequency offsets. 

The frequency offset compensator disclosed advantageously 
receives the estimate of the DC offset from the DC estimator and uses 
' the estimate of the DC offset to produce the estimate of the frequency 
offset, thus providing an improved estimate of the frequency offset. 

20 With reference to FIG. 2 one embodiment of a DC and frequency 

offset estimation and compensation system 200 comprises an input 202 
for receiving data packets 204 having a preamble 206. The input 202 is 
coupled to a DC offset estimator 210, a frequency offset estimator 212, 
and a compensator 214, and as the preamble 206 is received, it is 

25 provided concurrently to a DC offset estimator 210, a frequency offset 
estimator 212, and a compensator 214. The DC offset estimator 210 
produces an estimate of DC offset and is coupled to provide the 
estimate to both the frequency offset estimator 212 and the 
compensator 214. The frequency offset estimator 212 uses the estimate 

30 of DC offset and is coupled to provide an estimate of frequency offset to 
the compensator 214. The compensator 214 receives estimates of both 
the DC and frequency offsets, and is coupled to an output 216, which 



10 

provides a corresponding DC and frequency compensated data packet 
218. 

The system 200 performs concurrent estimation of the DC and 
frequency offsets when the preamble 206 is being received. With 
5 reference to FIG. 3A, this provides the opportunity for compensation to 
also be performed during reception of the preamble 206. Preamble 
reception time 305 is the time during which the preamble 206 of the 
data packet 204 is received. The preamble time 305 can be divided into 
a first portion, which will be referred to as estimation time 310, and a 

10 second portion, which will be referred to as compensation time 315. 
During the estimation time 310, the DC offset estimator 210 and the 
frequency offset estimator 212 operate concurrently, as shown by the 
blocks 320 and 332. Then, during the compensation time 315, the 
compensator 214 operates as shown by blocks 330 and 335 to 

15 compensate the received preamble 206 of the data packet 204 for both 
DC and frequency offsets, and can begin to provide the compensated 
data packet 2 18 by the end of the preamble reception time 305. 

Alternatively, with reference to FIG. 3B, the DC and frequency 
offset estimation can be performed during the preamble reception time 

20 305, and the DC and frequency offset compensation can be performed 
later i.e. after the preamble reception time 305, during a post 
processing compensation time 320. Of course this approach 
necessitates the use of a buffer to store the received data packet 204, 
and subsequently apply the post processing compensation to the stored 

25 data packet. 

With reference now to FIG. 4, the operation 400 of the system 200 
starts 405 with determining 410 when reception of the data packet 204 
has started, and when the determination is positive, providing 415 the 
received portion of the data packet 204 to the DC offset estimator 210, 

30 the frequency offset estimator 212, and the compensator 214. 

The DC offset estimator 210 estimates 420 DC offset of the data 
packet 204 from a part of the preamble 206, and provides 425 an 
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estimate of the DC offset to the frequency offset estimator 212. The 
frequency offset estimator 212 uses the estimate of DC offset to 
estimate 430 frequency offset of the data packet 204, and produces an 
estimate of the frequency offset 
5 The compensator 214 is provided 435 with the estimates of the DC 

offset and the frequency offset, and the compensator 214 uses the DC 
and frequency offsets to compensate 440 the data packet 204. Here, 
compensation is performed as the data packet is received, and the 
compensator 214 continues to compensate 440 the data packet 204 

10 until a determination 445 is made that the compensation of the data 
packet 204 is complete, and the compensated data packet 218 has been 
provided 450. The operation 400 then returns to determining 410 the 
start of reception of another data packet. 

With reference to FIG. 5 another embodiment of a DC and 

15 frequency offset estimation and compensation system 500 comprises an 
input 502 for receiving data packets 204 having a preamble 206. The 
input 502 is coupled to a DC offset estimator 505 and a DC offset 
compensator 507, The DC offset estimator 505 produces an estimate of 
DC offset from a portion of the preamble 206 as it is received, and has 

20 an output that is coupled to provide the estimate to the DC offset 
compensator 507. The DC offset compensator 507 uses the estimate of 
DC offset and the data packet 204 to produce a portion of a DC offset 
compensated preamble (not shown). 

The DC offset compensator 507 has an output which is coupled 

25 to provide the portion of the DC offset compensated preamble 
concurrently to a residual DC offset estimator 510 7 a frequency offset 
estimator 512, and a compensator 514. As the portion of the DC offset 
compensated preamble is received, it is provided concurrently to a 
residual DC offset estimator 510, a frequency offset estimator 512, and 

30 a compensator 514. The residual DC offset estimator 510 produces an 
estimate of residual DC offset and is coupled to provide the estimate to 
both the frequency offset estimator 512 and the compensator 514. The 
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frequency offset estimator 512 uses the estimate of residual DC offset, 
and is coupled, to provide an estimate of frequency offset to the 
compensator 514. The compensator 514 receives estimates of both the 
residual DC and frequency offsets, and is coupled to an output 516, 
5 which provides a corresponding residual DC and frequency 
compensated data packet 518. 

The system 500 performs sequential DC offset estimation and 
compensation followed by concurrent estimation of residual DC and 
frequency offsets when the preamble 206 is being received. With 
10 reference to FIG. 6A, this provides the opportunity for further offset 
compensation to be performed during reception of the preamble 206. As 
before, the preamble reception time 305 is the time during which the 
SJ preamble of the data packet 204 is received. The preamble time 305 can 

{U be divided into a first portion, which will be referred to as estimation 

15 time 605, a second portion, which will be referred to as DC offset 
compensation time 610, a third portion, which will be referred to as 
concurrent compensation time 615, and a fourth portion, which will be 
referred to as DC and frequency offset compensation time 610. During 
m the estimation time 605, only the DC offset estimator 505 operates as 

20 shown by the block 625, and during the DC offset compensation time 
610 only the DC offset compensator 507 operates as shown by block 
630. However, during the concurrent estimation time 615 both the 
residual DC offset estimator 510 and the frequency offset estimator 512 
operate, as shown by blocks 635 and 640. Finally, during the DC and 
25 frequency offset compensation time 620, the compensator 514 
compensates the portion of the DC offset compensated preamble, and 
can begin to provide the compensated data packet 518 by the end of the 
preamble reception time 305. 

Alternatively, with reference to FIG. 6B, the DC estimation and 
30 compensation, and the concurrent estimation of DC and frequency 
offsets can be performed during the preamble reception time 305, and 
the residual DC and frequency offset compensation can be performed 
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later i.e. after the preamble reception time 305, during a post 
processing compensation time 320. As before, this approach 
necessitates the use of a buffer to store the received data packet 204, 
and subsequently apply the post processing compensation to the stored 
5 data packet. 

With reference now to FIG. 7, the operation 700 of the system 500 
starts 705 with determining 710 when reception of the data packet 204 
has started, and when the determination is positive, providing 715 the 
received portion of the data packet 204 to the DC offset estimator 505 
10 and the DC offset compensator 507. The DC offset estimator 507 
determines 720 an estimate of DC offset of the data packet 204 from the 
y received portion of the data packet 204, and when a determination 725 

y3 is made that a predetermined estimation time 605 has elapsed, then the 

!*! DC offset estimator 507 provides 730 an estimate of the DC offset to the 

jy- 15 DC offset compensator 507. The DC offset compensator 507 uses the 
estimate of DC offset and the received portion of the data packet 204 to 
compensate the received portion of the data packet, and subsequently 
produce a DC offset compensated portion of the data packet 204 after a 
Ml determination 740 that the DC offset compensation time 610 has 

?y 20 elapsed. 

The DC offset compensated portion of the data packet 204 is then 
concurrently provided 745 to the residual DC offset estimator 510, the 
frequency offset estimator 512, and the compensator 514. The residual 
DC offset estimator 210 estimates 750 residual DC offset of the 
25 compensated portion of the data packet 204, and provides 755 an 
estimate of the residual DC offset to the frequency offset estimator 512. 
The frequency offset estimator 512 uses the estimate of the residual DC 
offset to estimate 760 frequency offset of the compensated portion of the 
data packet 204, and produces an estimate of the frequency offset. 
30 The compensator 514 is provided 765 with the estimates of the 

residual DC offset and the frequency offset, and the compensator 214 
uses the residual DC and frequency offsets to further compensate 770 
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the compensated portion of the data packet 204. Again, compensation 
is performed as the data packet 204 is received, and the compensator 
514 continues to compensate 770 the data packet 204 until a 
determination 445 is made that the further compensation of the data 
5 packet 204 is complete, and the compensated data packet 518 has been 
provided 780. The operation 700 then returns to determining 710 the 
start of reception of another data packet. 

The frequency offset estimator 212 and 512, in accordance with 
the present invention, will now be described based on a DCA receiver 
10 with OFDM modulation. 
l* As an illustrative aid, a model is provided with the following 

notation: 

f c is inter-carrier spacing of OFDM sub-carriers; 
f a is actual frequency offset of OFDM sub-carriers; 
} j? 15 € is the relative frequency offset i.e. the ratio of f a to f c \ 

i P is the number of sub-carriers used; and 

{Tj T= 1 / (P f c ) as the sampling interval. 

fl FIG. 8 shows the structure of the preamble of data packets that 

y ? 

O will be used in a simulation with the model, where the preamble 

iy 20 complies with the IEEE 802.11a standard. The preamble comprises M 

slots each having an inter-slot index m, and N intra-slots each having 
an intra-slot index n. The model will be used in a simulation, with P=64 
and N=16 on a communication channel 

A transmitted preamble can be defined as 
25 x(m,n) A \x(m,n)\exp(ja(m,n)) (1) 

The complex channel gain as: 

A A |A|expt/0) (2) 

And, as a result of the DCA, the DC offset as: 

8 A p|«p(/A (3) 
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The assumption that a time invariant complex gain A and DC 
offset 5 is valid, when the coherence time is relatively large as compared 
to the packet length, and when the gain controller at the receiver is held 
constant. 

The preamble is transmitted via a communication channel with a 
frequency offset defined as follows: 

r= 2 #« J 1 ™ 

(Pfc) (P) 

With reference to FIG. 9 the model 900 of the communication 
channel takes into account the various parameters of the 
communication channel that affect transmission of data packets. These 
include frequency offset 905 represented by exp(j^), complex gain 910 
represented by A = |A|exp(j0), direct current offset 915 represented by 
<5 = [5|exp(7jB), and noise 920 represented by v(m,n), which affect the 

transmitted preamble 925 represented by x(m,n) and results in a 
received preamble 930 represented by z(m 5 ft). Further, including 
additive white Gaussian noise (AWGN) in the communication channel 

with variance <7 V , the received preamble 930 is defined as follows: 

z(m 9 n) = y(m,n) + v(m,n) (4) 
y(m,ri) - Ax(m,ri) exipijyk) + d 

= \Ax(m y n)\exp[j(Yk + a(m,n) + 0] + |<^exp(/0) (5) 

PRIOR ART WITH NO DC OFFSET 
An article titled "Frequency Offset Synchronization and Channel 
Estimation for OFDM Based Transmission" by Song et al. in IEEE 
Communications Letters Vol. 4 No. 3 March 2000, discloses frequency 
offset estimation for a communication channel in the presence of 
AWGN, and where there is no DC offset. The article teaches that when 
the signal to noise ratio (SNR) is high, the frequency offset estimator is 
not biased. 
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This may be expressed, thus. 

E[e]-s = 0 

Then, the variance of the frequency offset estimator can be shown 
as follows: 

— 1 A <J 2 (M,N) (6) 

N2n) ( M -l) 2 N SNR= V ' W 

For unbiased estimators, the variance of the estimator is the mean 
square error (MSE) of the estimate. 

Another article titled "A Technique for Orthogonal Frequency 
Division Multiplexing Frequency Offset Correction" by Moose in IEEE 
q 10 Transactions on Communications Vol. 42 No. 10 October 1994, 
P discloses a maximum likelihood estimation algorithm. Using Moose's 

gg algorithm, the M repeated slots each with a length of N points is divided 

!., r, 

^ into two bigger slots, each of length MN/2, on the assumption that Mis 

fj : even. The variance of frequency offset estimate can then be obtained by 

JL 15 replacing M by 2 and N by MN/2 in equation (6), producing the 

yj following equation. 

| ct 2 (2,MV72)~ -=— - (7) 



yMNn) MN SNR/2 

Then assuming Mis even we have: 

\ M = 2 
>1, M = 4 
<1, M-6 5 8. 



cr 2 (2,^/2) ^ 8^-l 2 ) ^ 
<t 2 (M,N) M 3 



20 From the equation above, Moose's frequency offset estimator 

provides better performance than Song's in the way of lower variance for 
M =6, 8,.. etc. In addition, the improvement is marginal for smaller 
values of M. However, neither Moose nor Song considers a frequency 
offset estimator that takes into account the presence of a DC offset. 

25 
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PRIOR ART WITH DC OFFSET 
Now using the frequency offset estimator above with DC offset and 
equation (20), which will be derived later, the mean square error (MSE) 
of the frequency offset estimate, taking into account the DC offset, can 
be determined as follows; 

E{[s-E]} = E{[e-E(e) + E(E)-s] 2 } 

= E{[e-E(e)Y} + E{[E(e)-e]>} 

= E{[e-E(e)] 2 } + [E(e)-ef (8) 

This relates the MSE to the variance and bias of the estimator. 
When DC offset is present, the estimator is biased: 

m ~ * = ~ + ^ eXp0W ) ( 9 ) 

Thus, with the variance given in Equation (6), and if proper 
preamble is used, the MSE can be approximated as follows: 

(10) 

The accuracy of the result is validated by comparing equation (10) 
with the MSE of frequency offset estimates obtained from simulations 

for different values of 6 and \S[ using M = 4 at SNR = 20dB. With 
reference to FIG. 10, the MSE is normalized to £ and is consistent with 
the analytical expression. It can also be observed that the frequency 
offset estimate is sensitive to the presence of DC offset. 

DERIVATION OF EQUATION f20V 
Song's frequency offset estimator assumed the communication 
channel to be noisy and averages over a number of terms to improve the 
estimation. However, here an expression for the frequency offset 
estimation is first derived assuming that the channel is noiseless and 
no averaging is carried out. 



E[(e-ef] - a\M,N) +[^ arg(l + |«f ex 9 (jN Y j) 
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With reference to FIG. 11, which illustrates a frequency offset 
estimator based on r yy (m,m + l;ri), where r yy {m,m + \\n) is defined as 

the multiplication of discrete point y(m,n) with its conjugate N samples 
or the equivalent of one slot apart in time. Hence, 
Tyy (m,m + l;n) = y(m,n)y* (m + l,n) 

= |Ax(m,ra)| 2 exp(- jNy) + 1<5| 2 + 

\A8x(m,n)\{tx])(j£(m 9 n)) + exp(- j£(m + l,n))} 
= \Ax{m,nf exp(-//Vf) + \S\ 2 + 
\A8x(m t n)\&xp(j^(m 9 n)){l + exp[-j(£(m +l 9 n) + $(m 9 n))]} 9 (11) 
where £(m,n) A y(Nm + n)+a(m,n) + d-p. (12) 
Note that £(m + l,ri) + £(m,n) = 2%{m,ri) + Ny, and using the identity 
1 + expO'20) = 2cos(0)exp(/0), the last term of equation (1 1) can be 
written as follows: 

\A8x(m,n)\zxp(j£(m,n)){l + exp[-;(2|(m,n) + Ny)]} 
= \A5x(rn,n)\ • 2 cos(£(m,rc) + ATy/2) • exp(- jNy 12) (13) 

Now substituting equation (13) in equation (11) and taking out the 
factor exp(~jNy) the following equation is obtained. 

ryy (m,m + l;rc) = ^(-jNy)^Axi^m 9 n)f + 15| 2 cxp(jNy) + 

0(m,tt)exp(;tfy/2)] (14) 
wftere 0(m,n) A 2\ASx(m,n)\zos(£(m,n) + Ny/2). (15) 
Estimating the frequency offset using one point, when 
8 = 0,(j>(m,n) = 0 9 and the phase of r yy (m,m + l;n) is -Ny, assuming 

|Atyj<#, the frequency offset y under noiseless conditions can be 
obtained by dividing the phase by -AT. This is Song's frequency offset 
estimator with the assumption that no averaging is performed. 

When DC offset is present, <5*0, the frequency offset estimate 
then becomes: 
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f = ~ arg[r w (m, m + l;n)] 

1 

= y-— arg[ Ax(m,n)\ 2 + 1(5| 2 exp(jNy) + 0(m,n)exp(;iVy/2)] 



1 

= y arp 1 + 

N 



\Ax(m,n) 



^Ax(m,n) 



(16) 



The last two terms in the argument are undesirable and cause 
errors in the frequency offset estimate f . While the first undesirable 

term can be used for estimation if |<5| 2 is known, the second undesirable 
term is difficult to determine in practice. 

Based on the model with Song's frequency offset estimator, where 
here the communication channel is assumed to be noiseless, and now 
with averaging to isolate the undesirable term, an alternate expression 
for the frequency offset estimation is now derived. 

With reference to FIG. 12, which illustrates another frequency 
offset estimator based on averaging of r yy (m,m + l;ri), the time average 

of r yy (m,m + l;ra)is defined as R yy {m,m^r\\n). To simplify analysis, and 
for the substitution x(m,n) = x(0,n) to be valid in the summation, 
consideration is only accorded to the case when N(M -1) terms are 
used for the averaging over equation (15), thus: 
1 



iV(M-l) ffi f 0n t 0 

| N-l 



M-2N-1 

X %ryy(m,m + l\n) 



" B=0 



1 



M-2N-1 

X I>(m,«)exp(./iVy/2) 



]_N{M-l)^ 
= exp(-jNy)[l + |<f exp(/Wr) + W(M ,N, Y )] (17) 

where the undesirable term that is difficult to estimate, using equations 
(12) and (15) is given as: 
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2jA < 5|exp(;iV7/2)V^ 1 i , A m , , 

= ^(m-i) ^ XWM<™(r*0*»+»)+ 

a(O,n) + 0-0+Ary/2). (18) 
We normalized the power of the received preamble when DC offset 
and noise is absent as 1, i.e. 



l -l\Ax(<>4 2 =1 ( 19 ) 

iV n=0 



The implications are that \A\ is assumed to be known, and more 

importantly, that the values |<5| 2 and Y(M,iV,y) are normalized to the 
received preamble power. In the absence of DC offset, the frequency 
offset estimate, which is valid for \Ny\<7Z assuming that 

arg[l + 1<5| 2 expC/iVr + "¥(M,N,y)\ « Ny is, 
?=~aig[R yy (m i m + l;n)] 

= y- -^arg[l + |<5| 2 e X p(yiVy) + V{M,N,y)1 (20) 

As in equation (16), the last two terms in the argument cause 
errors in the frequency offset estimate even in the absence of noise. 
However, assuming |6j is known, and if ¥(M JV,y) « 1 , a more 
simplified expression can be obtained for f , the frequency offset 
estimate. 

VALIDITY OF W(M,N,y)«l 
Support for the assumption that Y(M,iV» « 1 now follows, with 
determining the conditions that need to be satisfied for the assumption 
to hold true when If/ A d-fi + Nyllz [0,2^) varies for different data 
packets. A method is proposed to design the preamble of the data 
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packets such that the assumption is valid for small values of y but for 
different y/ 7 since in practice this is an unknown value. Three 
conditions, Condition A, Condition B, and Condition C, are deduced 
which lead to increasing support as to the validity of 
^{M.N.y) « lwhen y is small. In particular it will be shown that the 
preamble that conforms to the standards of IEEE 802.11a satisfies the 
Condition A. 

For simplification, equation (18) is re-written as follows: 
V(M,N,yy) = 2\AS[ Y'(M,tf,y,¥f) (21) 
where 

Y(M,N,w) 4 **& N - ffj } %(0,n)\ YcostyNm + yn + g(0,n) + y) (22) 



m=0 



n-1 



Using the identity Xcos(Ax + a) = cosec— sin 



f \ 
fix 



cos 

for the summation over index m, equation (22) becomes, 

((M-^yN^ 

\ " 2 """" " " 

i N-l 

^(N.yy/) 4 — X|*(0,n)|cos(i« + a(0,n) + i/) 



+ <3 



w./'i-r at \ exp(;iVy/2) /jflA 
x ¥ , (M,N,Y,w) = — J -^ L - - — ^cosec — 
V 7 ^ M-1 \2) 



SUl 



tp' A y/+yN(M-2)/2e [0,2k) 



^"(iV.y.vO (23) 

(24) 
(25) 



A more useful expression of equation (24) for observing the 
conditions when ^"(Af^ty/Ocan be made small by some choice of 
15 x(0,n)is given as follows: 

N-l 



W"(N 9 y,\f) = -j-B t \x(Q,n)\ exp(;» exp ja(0,n)) expQY) 



=— 9? exp(yy)Xx(0,n)exp(j>) 

^ L n=0 



(26) 



Summarizing, the following equation results. 
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Y(M,2V,y,vr) = 



2) 



SIT 



MyN\\ 

Jn 




AM 



exp(;vO X^(0,n)exp(;» 

n=0 



(27) 



A, Bound and Statistical Analysis of ,y,yO 

For each preamble received, all the parameters and signal at a 
given (m,n)can be considered to be fixed. However, when the next data 
packet is received, we have a new set of parameters and signal. The 
bound of *F ,f (N»y>¥0 will now be calculated for any given 
y/e [0,2n) and ye [-nlN.nlN). Then assuming y/to be a random 

variable, the mean and variance of *¥ n (N,y,l[/) is analyzed. 

Using equation (26) and that in general for any value of 6 y 
|9t[exp(jfl)(a + jfr)]|=|cos(0)a - sm(d)b\=\a + jb$cos(d+atm(b/a))\ (28) 

the bound of can be expressed as follows: 

N-l 



E^(0,n)exp(y» 



/!=0 



Amax| Y"(N,Y,V)\ 



(29) 



Assuming that Iff is a random variable with uniform distribution, 

1 

with the probability density function (pdf) given as — for [0,27r), the 

2n 

15 mean of x F M (A^,y 5 l//)is, 



N 
J_J_ 
#2;r 



2tt 



jV-1 



L«=0 



= 0 (30) 
The variance can be obtained using equation (28) as follows: 
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var[ ¥"(V0]=J 0 2, 1 ^\¥f^-dV 



In 



N-l 



IN' 



£x(0,n)exp(y>w) 



n=0 



maxl T"(iV,y,vOl 



(31) 



Thus, the above relates the bound of | T'j and the variance of 
Y"with respect to \f/. Therefore, if max^| T"(A^,y,vO| can be made 
close to 0, the assumption that ¥(M,N,y,V^) « 1 is valid. This is true 
when M — >°o for any given S and y/, but not necessarily true when 
N-*oo, Typically, not many slots will be available for averaging, so a 
preamble is designed such that it has a small bound. This can be done 
by scrutinizing the function max^| Y" (N,y,V)\ • 



B. Condition A 

When yn«l for n = 0,...,N-l, then it can be written 
thatexp0» ~ l,« = 0,.„,iV-l. Then from equation (29), 

1 



max| Y M (JV,7)|« 



N 



N-l 

X*(o,«) 

n=0 



(32) 



(33) 



For ¥(M,AT,7,^) « 1, the Condition A is denned as follows: 

N-l 

£*(0,n) = 0. 

71=0 

C. Condition B 

Now assuming that N is even, an alternative expression for 
equation (29) is: 



maxl Y"(iV,y,¥0| = — 
¥ 1 V A N 



NI2-1 



E*(0,2«)exp(yy2n) + x (p,2n + 1) exp(jy(2n + 1)) 



N 



n=0 
N/2-1 



£exp(;y2rc)[*(0,2rc) + x(0,2n + 1) Qxp(jy)] 



H=0 



(34) 
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Assuming that where this is a less strict assumption than 

that used in equation (32), 

x(0,2rc) + *(0,2n + l)expt/y)« jt(0,2rc)-hx(0,2n + l), which is set to 0, 
Thus, the Condition B, which also satisfies the Condition A, is defined 
as follows: 

x(0,2n) = -x(0an + l\n = 0,- • -,^72 - 1 (35) 
If the Condition B is satisfied, 



nmx|T-(^r,¥0| = ^|(l-exp(iy))| 



AT/2-1 

X^xp( i /y2n)jc(0,2n) 

n=0 



(36) 



max| Y"(AT,n^| = l|x(0,0)| 



D. Condition C 

From equation (34) the procedure can be repeated N' = iog 2 N, 

where Wis a positive integer, times on equation (36) to obtain the 
following expression, 

nfl-exp^)) (37) 

71=0 

The bounds of equation (37) are obtained by subjecting the 

preamble to a set of conditions, which collectively is known as the 
Condition C, and defined as: 

x(0,2n) = -x(0,2n + l), n = 0,---iV 72-1 

x{0,4n) = -x(0,2(2n + 1)), n = 0,- • -iV/4 - 1 

*(0,8n) = -x(QA{2n + 1)), n = 0,- • • N/S - 1 

:=; (38) 

x(0,(N/2)n) = -x(0,(N/4)(2n + 1)), n = 0,- • • N /(N/2) - 1 
jc(0,0)=-jc(0,;V/2) 

To obtain the preamble which follows the Condition C, the 
following method is proposed: for N = 2, the preamble, which is 
represented as a row vector, is x 2 = *(0,0)[1,-1], for N = 4 the preamble 
is x 4 =[x2,-x 2 ]; and in general, the preamble is x^ = [^n/2» - %/2] • 
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Thus, the preamble is properly defined by specifying one of the 
preamble elements, say jc(0,0). For reference, for N =16, the preamble 

is x(0,0)«[l -1,-1,1 -1,1,1-1,-1,1,1,-1,1,-1-1,1]. For simulation x(0,0)is 

set to equal 1. It should be noted that, 

v2 



N'-l 



N'-ll 



+ • 



n(l-exp(yy2' l ))=n;y2' l -(y2 n ) 12 
with the lowest order term given as, 



N' 



Thus, assuming that y<^VN, the lowest order term decreases at 
a rate of N' and makes T M (N>y ? tfO very close to 0, regardless of the 
value of iff. The assumption V < appears to be less stringent than 

that required for the Condition B, which assumes that y «1. 

Hence, in view of the analysis above, the assumption that 
^(MfNrf) « 1, and can therefore be ignored in equation (20) is valid. 



FIRST EMBODIMENT OF A FREQUENCY OFFSET ESTIMATOR 
Based on the assumption that ^(M,iV,7) can be ignored, 
equation (20) can be simplified so that the frequency offset estimate f is 

given by, y = y - ^ arg[l + 15| 2 exp(;A7y)] (39) 

An approximated form of equation (39) is derived as follows: 

y = y--^arg[l+ |<5| 2 cos(ATy) + y|6f sin(tfy)] 

~y-^arg[l + |<f +y|<f tfy] 



=y — -tan 1 



Y N 



[<f iVy 

2 



1 + 



=y(l + |f)" 1 . 



(40) 
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In the first approximation it is assumed that Nf is small, so that 
cos(JVy) = 1 and sin(A/y) ~ Ny, and in the second approximation it is 

\di 2 Ny 2 
assumed that u — T «l le. Nv«l + M , such that 

1 + |<5| 2 

\3i 2 Ny 

j. Note that the approximation is exact when 



tan -1 



\$ 2 Ny 



1+|<5| 2 J 1 + |61 

|<5| 2 = 1 , which results in the estimation being more accurate in the 
region when |(5| 2 = 1, and the simulation results bear this out. 



In practice the frequency offset estimate f is estimated from the 
phase R^. Thus, the frequency offset estimate for this embodiment, 
which will be referred to as FOEI, is defined as follows: 

i*A-larg(/y(l + 0 2 ) (41) 

With reference to FIGS. 13 and 14 a frequency offset estimator 
1300 and its operation 1400 starts 1402 with determining 1405 
whether receipt of a data packet 1310 at an input 1305 has started. The 
data packet 1310 has a preamble 1315, which is represented by 
samples z(m 9 n). The input 1305 is coupled to a frequency offset 
estimator module 1320 and to a DC offset estimator 1325, both of 
which receive the preamble 1315. The DC offset estimator 1325 
processes the preamble to determine 1410 and produce an estimate of 

DC offset power, that is represented by |<5| . Subsequently, an adder 

1330 adds 1415 a predetermined constant to the DC offset power |<5| . 
In the frequency offset estimator module 1320, a memory module 1335 
stores 1420 a predetermined number of the samples, represented by AT, 
and the N samples are then processed by a complex conjugation 
module 1340, that performs 1425 complex conjugation on the N 
samples. The result from the complex conjugation module 1340 and the 
preamble z(m,n) 1315, are then multiplied 1430 together by a 
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multiplier 1345 to produce an output, represented by r K (w,m + l;«). An 
averaging module 1350 averages 1435 r^(m,m + l;rc) over a 
predetermined number of samples, K, and the averaging module 1350 
provides an output represented by /f^(/n,/n + l;n). An argument divider 
1355 then determines the argument of R^m.m + l\n), and divides 1440 
it by -N to produce a resulting quotient. The output of the argument 
divider 1350 i.e. the resulting quotient, and the output of the adder 
1330 are then multiplied 1445 together by multiplier 1360, and the 
multiplier 1360 then provides 1450 a frequency offset estimate 7 to an 
output 1365. The operation 1400 then returns to determining receipt of 
another data packet, and repeats as described. 

SECOND EMBODIMENT OF A FREQUENCY OFFSET ESTIMATOR 
Again assuming that *F(Af can be ignored, equation (17) can 
be rewritten, as follows: 

Ryy (m 9 m + l;n) * exp(-jNy) + |<f (42) 
This is estimated from R u {m,m + \\ri), thus where implementation 
allows |<5| 2 to be removed from first before calculating the phase, 

then the frequency offset for this embodiment, which will be referred to 
as FOEII, is defined as follows: 

i?/A-larg(^-^) (43) 

With reference to FIGS. 15 and 16 a frequency offset estimator 
1500 and its operation 1600 starts 1602 with determining 1605 
whether receipt of a data packet 1510 at an input 1505 has started. The 
data packet 1510 has a preamble 1515, which is represented by 
samples z(m,ri). The input 1505 is coupled to a frequency offset 
estimator module 1520 and to a DC offset estimator 1525, both of 
which receive the preamble 1515. The DC offset estimator 1525 
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processes the preamble to determine 1610 and produce an estimate of 
DC offset power, that is represented by |<5| 2 . In the frequency offset 
estimator module 1520, a memory module 1535 stores 1620 a 
predetermined number of the samples, represented by N, and the N 
samples are then processed by a complex conjugation module 1540, 
that performs 1625 complex conjugation on the N samples. 

The result from the complex conjugation module 1540 and the 
preamble z(m,n) 1515, are then multiplied 1630 together by a 
multiplier 1545 to produce an output, represented by r zz (m,m + l;n). An 
averaging module 1550 averages 1635 r u (m,m + l;n) over a 
predetermined number of samples, K, and the averaging module 1550 
provides an output represented by R a (m,m + l;n). 

A subtracter 1555 combines the output fl a (m,m + l;n) from the 
averaging module 1550 and the output |of from the DC offset estimator 
1525, and an argument divider module 1560 then determines the 
argument of the result from the subtracter 1555, and divides 1645 the 
argument by -N to produce the frequency offset estimate y. The 
argument divider module 1565 then provides 1650 the frequency offset 
estimate y to an output 1565. 

The present invention, as described in the embodiments herein, 
advantageously provide a frequency offset estimator that takes into 
account an estimate of DC offset, and produces an improved frequency 
offset estimate. 

Using computer simulations with the number of discrete points in 
one preamble slot W = 16, for different SNR and 6, MSE is calculated 
and plotted for 1000 simulations against increasing DC offset power 
|5| 2 , to produce a series of graphs. 

With reference to FIGS. 17A-F the graphs show resultant MSE for 
both FOEI and FOEII across increasing DC offset power |<5| 2 , for the 
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three preamble conditions i.e. Condition A, Condition B and Condition 
C, when two preamble slots are used i.e. M=2. The relative frequency 
offset is 6 = 0.1 for FIGS. 17A and 17B; is £ = 0.5 for FIGS. 17C and 
17D; and is 6 = 1 for FIGS. 17E and 17F. The SNR for the graphs in 
FIGS. 17A, 17C and 17E is 20dB and can be compared with stronger 
signal conditions in FIGS. 17B,17D and 17F, where the SNR is 40dB. 
Similarly, with reference to FIGS. 18A-F the graphs show resultant 

MSE for both FOEI and FOEII across increasing DC offset power |<5| 2 , for 
the three preamble conditions i.e. Condition A, Condition B and 
Condition C, when four preamble slots are used i.e. M=4. The relative 
frequency offset is £ = 0.1 for FIGS. 18A and 18B; is £ = 0.5 for FIGS. 
18C and 18D; and is £ = 1 for FIGS. 18E and 18F. The SNR for the 
graphs in FIGS. 18A, 18C and 18E is 20dB and can be compared with 
stronger signal conditions in FIGS. 18B,18D and 18F, where the SNR is 
40dB. 

The FIGS.17A-F and 18A-F, show that preambles that satisfy the 
Condition A, Condition B and Condition C advantageously result in 
decreasing MSE for both FOEI and FOEII. The decrease is more 
significant for FOEII than FOEI, particularly at higher SNR values. 

With reference to FIGS. 19A-19C, the graphs show plots for a prior 
art FOE, referred to in the graphs as original FOE, and FOEII in 
accordance with the present invention, as described. Again using 
computer simulations with the number of discrete points in one 
preamble slot AT = 16 at SNR=20 dB for different £, MSE is calculated 
and plotted for 1000 simulations against increasing DC offset power 

J5| 2 , to produce a series of graphs for M=2, 4 and 6. 

The FOE of the present invention as described, advantageously 
provides substantially improved performance relative to the prior art 
FOE, particularly in the increasing presence of DC offset. 

This is accomplished by a frequency offset estimator that takes 
into account DC offset in determining a frequency offset estimate. The 
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frequency offset estimator can then be used with a DC offset estimator, 
that provides the DC offset estimate, and with a compensator, to 
provide a frequency offset estimation and compensation system. 

The frequency offset compensator advantageously receives the 
5 estimate of the DC offset from the DC estimator and uses the estimate 
of the DC offset to produce the estimate of the frequency offset, thus 
providing an improved estimate of the frequency offset. 

As the DC offset estimator and the frequency offset estimator 
operate concurrently to produce estimates of DC and frequency offsets 
10 from the preamble of received data packets, the concurrent operation 
h advantageously allows more time for both the DC offset estimator and 

the frequency offset estimator to produce the estimates of DC and 
frequency offsets. 

Thus, the present invention, as described provides a method and 
fp 15 apparatus for frequency offset estimation, and system utilizing same, 

t which overcomes or at least reduces the abovementioned problems of 

P 

y the prior art. 

I J. 

It will be appreciated that although only particular embodiments of 
the invention have been described in detail, various modifications and 
20 improvements can be made by a person skilled in the art without 
departing from the scope of the present invention. 
We claim: 
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